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Abstract 

The time evolution of a closed system of mean fields and fluctuations is Hamilto- 
nian, with the canonical variables parameterizing the general time-dependent Gaus- 
sian density matrix of the system. Yet, the evolution manifests both quantum de- 
coherence and apparent irreversibility of energy flow from the coherent mean fields 
to fluctuating quantum modes. Using scalar QED as an example we show how this 
collisionless damping and decoherence may be understood as the result of dephasing 
of the rapidly varying fluctuations and particle production in the time varying mean 
field. 



Mean field methods have a long history in such diverse atomic physics 

(Born-Oppenheimer), nuclear physics (Hartree-Fock), condensed matter (BCS) and 
statistical physics (Landau-Ginzburg), quantum optics (coherent/squeezed states), 
and semiclassical gravity. Because no higher than second moments of the fluctua- 
tions are incorporated, the mean field approximation is related to a Gaussian vari- 
ational ansatz for the wave function of the system. The broad applicability of the 
approximation, as well as the variety of different approaches to it in the literature 
makes it worthwhile to exhibit its general features unobscured by the particulars of 
specific applications. Accordingly, our first purpose in this Letter is to demonstrate 
the equivalence of the time-dependent mean field approximation to the general Gaus- 
sian ansatz for the mixed state density matrix p, and to underline its Hamiltonian 
structure. 

The Hamiltonian nature of the evolution makes it clear from the outset that the 
mean field approximation does not introduce dissipation or time irreversibility at a 
fundamental level. Nevertheless, typical evolutions seemingly manifest an arrow of 
time, in the sense that energy flows from the mean field to the fluctuations without 
returning over times of physical interest [Fig. 1]. Closely connected to this effective 
dissipation is the phenomenon of quantum decoherence flj, i.e., the suppression with 
time of the overlap between wave functions corresponding to two different mean field 
evolutions [Fig. 2]. This decoherence is the reason why quantum superpositions of 
different mean field states are difficult to observe in nature, and is crucial to un- 
derstanding the quantum to classical transition in macroscopic systems. Our second 
aim in this Letter is to present an explicit example of a quantum field theory (scalar 
QED) treated in mean field approximation where these effects are observed, and to 
provide a clear physical explanation of the behavior in terms of dephasing of the 
fluctuations, i.e., the averaging to zero of their rapidly varying phases on time scales 
short compared to the collective motion of the mean field(s). 

To expose the general structure of the time-dependent mean field (or Gaussian) 
approximation consider first a one-dimensional harmonic oscillator with Hamiltonian, 



where the frequency u(t) is a smooth function of time, otherwise unspecified for the 
moment. The most general Gaussian ansatz for the mixed state normalized density 
matrix may be presented as 
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in the coordinate representation. The five parameters (q,p, £, rj, p) of this Gaussian 
may be identified with the two mean values, q = (q) = Tr(gp), p = (p) = Tr (pp), and 
the three symmetrized variances via 

((<? - qf) = £ 2 > ((pq + qp- W)} = 2£v , 

(ip-p) 2 ) = V 2 + ^r- (3) 



The one anti-symmetrized variance is fixed by the commutation relation, [q, p] = ih. 
The parameter p measures the degree to which the state is mixed: Tr p 2 = pr l < 1, 
the equality holding for pure states. If the state is pure, p = and only two of 

the three symmetrized variances in (^j are independent. 

The Gaussian ansatz for the density matrix is preserved under time evolution. In 
the Schrodinger picture p evolves according to the Liouville equation, p = —i[H,p]. 
Substitution of the Gaussian form (§) into this equation with Hamiltonian (jl]) and 
equating coefficients of x, x', x 2 , x' 2 and xx' gives five evolution equations for the five 
parameters specifying the Gaussian, 

q = p; p=-u 2 (t)q 

and fx = 0. Since p, is a constant and the von Neumann entropy, — Tr plnp of the 
state (0) is a (monotonic) function of p alone, this quantity is also a constant of 
the motion. Evolution equations for the diverse applications of the time-dependent 
mean field approximation reduce to (multiple copies of) equations of precisely the 
general form of (£|) , with uo 2 (t) a different self-consistently determined function of the 
coordinates and time, depending on the application. This establishes the equivalence 
between mean field methods and Gaussian density matrices for all evolutions of the 
form of Eqns. @. We give explicit examples below. 

An essential property of the evolution equations @ is that they are Hamilton's 
equations (hence, time reversible) for an effective classical Hamiltonian M, with rj 
playing the role of the momentum conjugate to £, 

H eff (q, p; e, i) = Tr(p# ) = \{f + if) + V eff , (5) 

and V e ff(q, £) depending on the particular form of u 2 (q(t), £(t); t). For example, if 
the original system is an anharmonic double well with quantum Hamiltonian 

H{q,p) = \p 2 + \{q 2 -v 2 ) 2 , (6) 

then the (large N) mean field equations of motion are identical to Eqs. @) with 
u 2 {t) = X(q 2 (t) + £ 2 (t) - v 2 ). In this case, 



The resulting evolution equations (|]) are now quite non-linear (and chaotic). The 
last "centrifugal" term in the effective potential is a manifestation of the quantum 
uncertainty principle which prevents the Gaussian width £ from shrinking to zero. 
The unitary operator U(t) which effects the time evolution of the density matrix 

(1), 

p(t) = U(t)p(0)W(t) , U(t) = exp (-iJ o Hdt^j (8) 
is given explicitly in the coordinate basis by 

(x'\U(t)\x) = (2mhv(t))-* exp | (u(t)x 2 + v(t)x' 2 - 2xx') j (9) 

in terms of the two linearly independent solutions to the classical evolution equation, 
( d 2 A / u \ u(0) =v(0) = 1 

U + ^')(J^ 0i iU(o) = o. (10) 

The Gaussian dynamics may be expressed as well by means of a Fock representa- 
tion of the time dependent Heisenberg operators, 

q(t) = rf(t)q(0)U{t) = q(t)+af(t) + a!r(t) 

p(t) = U\t) P (0)U(t)=p(t) + af(t) + <Jr(t) (11) 

where [a, a*] = 1. The complex mode functions / satisfy the evolution equation ( [TU] ) 
and the Wronskian condition, 

//• - ./T '7, (12) 

showing that Gaussian time evolution is essentially classical, with h appearing only 
in the time independent condition ([12]) enforcing the quantum uncertainty relation. 
There is considerable latitude to redefine / by the Bogoliubov transformation / — ► 
cosh 7 e^ e+ ^/ + sinh7 e^ e ~^f* without affecting the Wronskian condition ([12]). Such 
Bogoliubov transformations form a noncompact Lie group, the metaplectic group 
Mp(2), which is a double covering of the symplectic group of classical Hamiltonian 
dynamics, Sp(2) = SU(1, 1), whose Lie algebra is generated by the three symmetric 
bilinears aa, a^a^ , and aa^ + a^a. If the a and operators are appended to these 
three, the algebra again closes upon itself, forming a five parameter group IMp(2) 
M. The unitary evolution (|8|), (|9]) of the Gaussian density matrix (0) is an explicit 
representation of this group. 

The group structure can be exploited to choose a basis in which all expectation 
values vanish, except 

(a j a) = (aa j ) - 1 = JV > . (13) 
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The Gaussian density matrix is diagonal in the corresponding a)a time-independent 
number basis, 

(n \p\n) = —— —— , 14 

n + 1 + V 

with /j = 2iV + 1 and £ 2 (t) = p\f(t)\ 2 . Upon identifying // = cot\i{huj/2kT), the 
diagonal form (K^f will be recognized as a thermal density matrix at temperature T. 
The pure state Gaussian wave function (/i = 1) corresponds therefore to a coherent, 
squeezed zero temperature vacuum state. The smoothness of the finite temperature 
classical limit Hp, — > 2kT/u as h — > 0, p — > oo shows that quantum and thermal 
fluctuations are treated by the mean field approximation in a unified way. 

By making another IMp(2) group transformation it is always possible to diago- 
nalize ([]]) at any given time, bringing the quadratic Hamiltonian into the standard 
harmonic oscillator form, H osc = ^ (acfi + a)a*j with a time dependent. This time 
dependent basis is defined by the relations, 

q(t) = af + a* f* p(t) = —itoaf + iua/f* 

/(t > = Sy^Hi'^ ) ■ (15) 

in place of (PH). In the a) a number basis, p is no longer diagonal, (a), (da), etc., are 
non-vanishing, and N = (a) a) ^ N in general, becoming equal only in the static case 
of constant uo. 

If u(t) varies slowly in time, an adiabatic invariant may be constructed from the 
Hamilton- Jacobi equation corresponding to the effective classical Hamiltonian fl5|). 
By a simple quadrature we find the adiabatic invariant, 

W (H) fi 
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Since N is time independent, N(t) is an adiabatic invariant of the evolution. On the 
other hand, the phase angle conjugate to the action variable W varies rapidly in time. 
Since the diagonal matrix elements of p in the iV basis are independent of this phase 
angle, they are slowly varying, whereas the off-diagonal matrix elements of p in this 
basis (which depend on the phase angle) are rapidly varying functions of time. If we 
are interested only in the effects of the fluctuations on the more slowly varying mean 
fields it is natural to define an effective density matrix p e ff{t) by time-averaging the 
density matrix @, thereby truncating p to its diagonal elements only, in the adiabatic 
iV basis Q. Clearly, for this truncation to be justified there must be very efficient 
phase cancellation, i.e. dephasing, either by averaging the fluctuations over time or 
by summing over many independent fluctuating degrees of freedom at a fixed time. 
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Obtaining the general form of the diagonal matrix elements of p in the N basis is 
straightforward but the result is rather unwieldy. Here we restrict our attention to 
the case of a pure state with vanishing q mean field. Using the methods of Ref. [Q, 
one finds simply 



(n = 2t\p\n = 21) 



, ' 2 l i\ 
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9=p=0 



-sech7 tanh 7 (17) 



with pnn = for n odd and j(t) the parameter of the Bogoliubov transformation 
between the a and a bases. It is given explicitly by sinh 2 7 = N = \f + iujf\ 2 /2hui. 

Decoherence is also addressable within the same framework. Consider the case 
where u(t) is a function of one external mean field degree of freedom A(t). If only 
the evolution of A is of interest, then the fluctuating modes described by f(t) may 
be treated as the "environment." To solve for the evolution of the reduced density 
matrix of A, one needs to compute the influence functional. This is a functional of 
two trajectories Ai(t) and A 2 (t) (corresponding to two different evolution operators 
Ui(t) and Uv(t)), and is defined by 

F 12 (t) = exp(ir 12 (t)) = Tr ([/i(t)p(0) U%(tf 



Explicit evaluation may be carried out using (|Sp. Restricting again to pure states 
with vanishing q mean fields we find 

ih 
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(19) 



in terms of the two sets of mode functions f\(t) and f 2 (t) which satisfy (^0|) and 
(P^). This T 12 is precisely the closed time path (CTP) effective action functional 
which generates the connected real time n-point vertices in the quantum theory 
For a pure initial state, the absolute value of measures the overlap of the two 
different evolutions at some time t, beginning with the same initial \tp(0)). In mean 
field theory, instead of evaluating T 12 for two arbitrary trajectories, the evaluation 
is over trajectories determined by the self- consistent evolution of the closed system, 
beginning with two different initial mean fields. 

For an explicit field theoretic example consider scalar QED with no scalar self- 
coupling. In the large N limit, the evolution of electric fields and charged matter field 
fluctuations may be described in the self-consistent mean field or Gaussian approxi- 
mation pi . For a spatially homogeneous electric field in the z direction, in the gauge 
A = A(t)z, the time evolution equations read 

Mt) = (j z (t)) = I E (*- - eA(t)) \ft(t)\\2N(k) + 1) 
a- 
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fi(t) = (20) 



with the Wronskian condition ([12]) holding for every discrete wave number k in the 
finite volume V. In field theory there are an infinite number of quantum fluctuating 
modes ft of the charged scalar field, each varying rapidly in time with its own char- 
acteristic frequency, = (k — eA(t)) 2 + m 2 . The Gaussian density matrix is an 
infinite product of Gaussians each of the form (|2p and there are an infinite number 
of Bogoliubov parameters ^(k;t), one for each k. The equations of motion (p0|) are 
again Hamiltonian in structure. In this case the effective classical Hamiltonian, 
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+ ^E[l4l 2 + ^II4l 2 ](2iv(fc) + i) (21) 



describes charged particle production in the electric field E = —A by the Schwinger 
mechanism and the effects of the current {j z (t)) generated by these charged particles 
back on the electric field, through the semiclassical Maxwell equation in (|20D . The 
mean value of the scalar field itself is zero so that we may use the expressions ( |TTD 
and (|19[) for the effective density matrix and decoherence functional of the charged 
field fluctuations. The values, p 2 e given by ( |17D are then the probabilities of observing 
I charged particle pairs in the adiabatic N basis. The diagonal matrix elements of p 
for odd n vanish because particles can only be created in pairs from the vacuum. 

In this specific model we present numerical results (in 1 + 1 dimensions with 
vacuum initial conditions) on damping and decoherence of the mean electric field in 
Figs. 1 and 2. The plasma oscillations of the mean electric field on timescales long 
compared to cu^ 1 are clearly seen, as well as the apparently irreversible flow of energy 
from the electric field towards the charged particle modes. In Fig. 1 we plot both the 
Boltzmann entropy and the von Neumann entropy of the effective density matrix, for 
comparison. These are defined by 

S B = £{(jV(£) + l)ln(iV(jfe) + l) - N(k) In N(k)} 

k 

oo 

S eff = -TtPeff In Peff = - X) £ pu(k) In p 2e (k) (22) 

k ^=o 

respectively. Both display general increase during intervals of particle creation [0, 
when the electric field is sufficiently strong for the Schwinger pair creation mechanism 
to be effective. Neither quantity is a strictly monotonic function of time (no H- 
theorem). Since the charged particle modes ft interact with the mean electric field but 
not directly with each other, the effective damping observed is certainly collisionless, 
and may be understood as due to dephasing similar to that responsible for Landau 
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Figure 2: Absolute value of the decoherence functional as a function of time. The two 
field values are E and E — A. The top figure shows (for fixed A) the sharp dependence of 
decoherence on particle production when \E\ > .2E C . The second illustrates the relatively 
milder dependence on A. 
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damping of collective modes in classical electromagnetic plasmas. The entropy S e ff 
of the effective density matrix provides a precise measure of the information lost 
by treating the phases as random. The Boltzmann "entropy" would be expected 
to equal S e ff only in true thermodynamic equilibrium, which is not achieved in the 
collisionless approximation of Eqns. (|20|) . Otherwise we see from Fig. 1 that Sb 
generally overestimates the amount of information lost by phase averaging. 

That decoherence is closely related to the same dephasing of the particle modes is 
seen most clearly by comparing the absolute value of F 12 for different initial electric 
fields. Decoherence is very slow for electric fields less than the Schwinger pair produc- 
tion threshold but becomes very rapid above it ||]. This shows the strong dependence 
of the decoherence process on the same particle production by the mean field. 

The authors wish to acknowledge several helpful discussions with F. Cooper and 
A. Kovner. 
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